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ABSTRACT 


This report describes a computer program which calculates the viscous 
effects on the lift and pressure distributions of arbitrary three-dimensional 
bodies. The program is a combination of a panel method which computes the 
potential flow about arbitrary three-dimensional lifting configurations, and 
a two-dimensional boundary-layer method, which calculates the viscous-effects. 
These effects are applied to the three-dimensional shape in a "strip-theory" 
sense and the resultant viscous lift and pressure distributions are produced. 
Two methods of simulating the boundary layer in the final potential-flow 
calculation are available: 1) addition of the displacement thickness to the 
original shape, and 2) defining a surface blowing distribution on the body. 


The computer program is written in Fortran IV for the IBM 370 systems. 
16 temporary external units are used for storage. The region size needed to 
execute the program is about 360K bytes, but this is a direct function of 
the number of elements defining the configuration. 


Also presented in this reportis a detailed description of the program 


logic, complete instructions for executing the program, and a sample case. 
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1.0 INTRODUCTION 


A procedure has been described [4] for the calculation of the viscous effects 
on the lift and pressure distributions of arbitrary three-dimensional configu- 
rations. This report is the user’s manual for the computer program for this pro- 
cedure which combines Hess' three-dimensional low-speed panel method [1] with 
Cebeci's two-dimensional boundary-layer method [2],[3]. The program automatic- 
ally corrects the calculated potential-flow field for the displacement effect of 
the boundary layer. The corrections are made along streamwise strips of the 
lifting configuration where the boundary layer is assumed to be two-dimensional. 


The logic of the program can conveniently be described in three steps. In 
step 1, the geometric quantities which define the configuration and the desired 
freestream flow conditions are input,and velocities and pressures are calculated 
at the control points of the body. The second step then determines the stagna- 
tion point of each lifting strip and, together with the previously calculated 
surface velocities, input quantities are generated for use in the boundary-layer 
program where a distribution of displacement thickness is determined. In the 
final step the calculated displacement thickness is used to simulate the boundary- 
layer effect. This effect may be simulated by either the addition of the dis- 
placement thickness to the original body coordinates to generate a new shape 
(referred to as the surface displacement method) or by the generation of a new 
set of source densities (blowing method). In either case a final potential-flow 
calculation which includes the viscosity effects is made. 


The program has been prepared for use on the IBM 370 series digital com- 
puters and is written in Fortran IV. It makes use of overlays and requires 16 
external storage units. The core size requirements depend on the number of ele- 
ments that define the configuration but about 360 bytes may be needed. Single 
precision arithmetic is used. The time for a computer run is a function of 


several variables (element number, lifting strips, boundary-layer simulation 
methods, etc.) but a "typical" case of 500 lifting elements should take about 
15 minutes. 


The purpose of this report is to outline the use of the program and, in 
particular, to provide complete input instructions and a detailed explanation 
of the output. The technical background of the method for which this program 
is written and a discussion of some calculated results is presented in reference 4. 


2.0 PROGRAM DESCRIPTION 


This program couples a three-dimensional lifting potential-flow program 
{1] and a simplified version of a finite-difference program for calculating 
compressible laminar and turbulent boundary layers [2], [3] with a program 
which simulates the effect of the boundary layer on the flow. Two alternate 
subroutines allow the boundary layer to influence the potential flow by a dis- 
placement of the surface coordinates or by a simulated surface blowing. The 
program first computes an initial potential-flow solution, then solves two- 
dimensional boundary-layer equations in streamwise strips, to generate a dis- 
tribution of displacement thickness over the body surface. This displacement 
thickness is used either to generate a new set of geometry coordinates or to 
calculate a simulated blowing velocity over the original body surface. The 
final potential-flow calculation is then made with the chosen houndary layer 
procedure producing the simulated viscous flow. 


2.1 Potential-Flow Calculation 


This part of the program is similar to that described in ref. [1]. It 
takes the body geometry coordinates, forms the velocity matrix, solves for the 
source density distribution, and computes the pressure coefficient at each of 
the element control points. 


2.1.1 Subroutines and Their Functions 


Following is a summary of the routines and their functions: 


Name of 

Routines Functions 

BDFORM Directs the input to form the body geometry. 
INPUT Accepts initial card input of the body points. 
DINPUT Takes in the body points with 6* added on. 
LIFT Forms the lifting elements. 

NOLIFT Forms the nonlifting elements. 

VFORM Generates the velocity matrix. 

VFMNLF Forms the nonlifting velocities. 


Name of 
| Routines Functions 
VFMLFT Forms the lifting velocities. 
AFORM Generates the normal velocity matrix to be solved for density 
distribution. 
MATSOL Directs the program to use initial matrix solution, or to solve 
more right sides. 
COLSOL Solves the columnwise matrix. 
SOLMOR Solves further right sides. 
COMFLO Computes the final combined flows. 
PKUTTA Solves the bound vorticity coefficients for the equal-pressure 
Kutta conditions. 
PRINT Generates the final output. 
HEADER Generates header label. 
LDAY 


2.1.2 Program Notes 


The body geometry input consists of sections and strips. Each section is 
divided into a number of strips, and each strip contains a certain number of 
elements (on-body elements and wake elements). These elements are defined by 


their corner points which are input along adjacent N-lines, two of which define 
a strip. A number of these N-lines compose a section as illustrated below: 


A SECTION 
STRIPS 


LIFTING STRIP 
OF ELEMENTS 


ON BODY 
ELEMENTS 


CORNER POINTS 


TRAILING EDGE | ELEMENTS THAT DEFINE 
wr ELEMENTS 
ELEMENTS 


Details of the input are available in reference 5. 


For the lifting section of the body, the elements are formed into strips 
which are aligned in a direction along the streamflow. Each strip must also 
have elements off the body surface which will define the trailing vortex wake. 
Calculation of the flow about nonlifting bodies, e.g. a fuselage, may also be 
calculated in a similar manner. In this case, however, no wake representation 
is necessary, and the elements need not be in streamwise strips. Further, no 
boundary-layer displacement corrections are made on the nonlifting sections. 
Calculations for configurations that have both lifting and nonlifting sections 
may be made in one computer run, i.e. a wing on a fuselage. 


Other special features of the program include the use of "extra strips" 
to carry wing vorticity through the fuselage and "ignored elements" to handle 
such problems as pylon-wing intersections. Section 6.8 of reference [1] pro- 
vides further information on this procedures. 


2.2 Two-Dimensional Boundary-Layer Calculation 


2.2.1 Program Notes 


A modified version of the Douglas finite-difference method for calculat- 
ing compressible laminar and turbulent boundary layers [2], [3] is used here 
with a simplified input subroutine. The only inputs required are: (1) ref- 
erence velocity in ft/sec, (2) Reynolds number/ft, and (3) the transition 
locations. The reference velocity and Reynolds number are input only once 
for the whole case, but transition locations have to be specified for the 
lower and upper surfaces separately on each strip. The user has the option 
of inputting the desired transition location number, or he may let the pro- 
gram calculate a transition location. Thus, an eight strip wing case requires 
8 transition locations to be specified for the lower surface, and 8 transition 
locations for the upper. 


The velocity and arc-length input needed for the boundary-layer program is 
determined inside the program. These parameters are stored on external units in 


the simulation program and passed to the boundary-layer program, 


In order that the coordinate and velocity data from the potential flow 
part of the program is compatible with the input and output of the boundary- 
layer part of the program, a constant, FC (foot converter) must be input 
to the potential-flow program. The boundary-layer program is dimensional and 
calculations are done in terms of feet, and the displacement thickness 6* is 
output in feet. Reference velocity and Reynolds number are also in ft/sec 
and ft"! respectively, as stated before. 


Two examples may help to clarify the situation. A straight untapered 
wing is given with the chord equal to 1.0 foot and a test Reynolds number 
equal to 1.6 x 10°. The coordinates are input as they stand with FC = 1.0 
and Re = 1.6 x 10°. If the reference velocity is not known, any small value, 
j.e. 100 ft/sec, may be chosen. Since the program calculates only incompres- 
sible flow, the value of the reference velocity does not affect the results. 


Now assume that a swept tapered wing is defined with the root chord equal 
to 42 inches and the tip chord equal to 21 inches. The Reynolds number is 
18 x 10° based on the mean aerodynamic chord. First, the data must be con- 
verted into feet; this can be accomplished by inputting the data as given (in 
inches), and inputting FC = 12.0. In this case the mean aerodynamic chord 
is used as the average chord and the Reynolds number per foot is calculated 
based on this value. The mean aerodynamic chord is 32 inches which equals 
2.67 ft, thus the input Reynolds number/ft would be 18 x 10°/2.67 ft = 6.74 
Xx 10° +". Again the reference velocity is given or some arbitrary smal] 
value is used. 


2.2.2 Subroutines and Their Functions 


Routines Functions 

BOUNDL Controls program responsible for calling other boundary-layer- 
related subprograms. 

INPT Handles all input data to the boundary-layer program. 

EINF Calculates the transformed J-grid points. 

IVPF Generates initial velocity profile. 

FLPR Computes the fluid properties. 

EDVS Computes the eddy viscosity. 


Routines Functions 

SHFT Provides the initial guesses for each station. 

MOMX Finds the solution of the x-momentum equation. 

TRNS Computes the location of boundary-layer transition. 
SLOP Computes the derivative DYDX from X,Y input. 


OTPT Handles the results of the boundary-layer calculation. 


2.3 Simulation Programs 


2.3.1 Program Notes 

This nart of the program manipulates the data in and out of the potential- 
flow and boundary-layer programs. It first determines the location of the 
stagnation point of a lifting strip, then separates the lifting strip into 
two "boundary-layer" strips, starting at the location of the stagnation point 
and proceeding toward the trailing edge cr both upper and lower surfaces. It 
also organizes the velocity and arc length daca for each "boundary-layer" 
strip. Experience has shown that the boundary-layer program requires at least 
a 50-station solution for good accuracy. Accordingly, the program interpolates 
to obtain input at 50 prescribed percent arc length stations before entering 
the boundary-layer program. 


Upon returning from the boundary-layer program, the calculated displace- 
ment thickness distribution of each boundary-layer strip is back-interpolated 
to the original control-point locations. The upper surface and the lower 
surface boundary-layer data is then transformed to properly represent the 
displacement thickness associated with each lifting element. The program now 
branches to either the blowing method of boundary-layer simulation or the 
surface displacement method, depending upon the user's choice. 


If the blowing method is used, a new onset flow with the boundary-layer 
influence is generated. The program proceeds to solve for a new set of source 
densities, and the potential-flow solution with viscosity effects is obtained. 


If the surface displacement method is selected, a new set of input body 
coordinates which have the displacement thickness added to them is generated 
and fed back to the potential-flow program. Another potential-flow calculation 
is them made with the new input body coordinates. 
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In either case, a complete set of the new potential flow output js 
printed. 


2.3.2 Subroutines and Their Functions 


Routines Functions 


SETUP Links the potential-flow program with the boundary-layer 
program. 


BSETUP Processes the input and output of the boundary-layer program 
for blowing method. 


DSETUP Same as BSETUP except it is for the surface displacement method. 


CALCRS Calculates a new flow to be solved in potential-flow program 
(blowing method). 


MODXYZ Modifies the original x,y,z coordinates with 6* added on 
and generates new X.Y,Z data to be used in the final potential 
flow run (surface displacement method). 


Handles the two-dimensional data interpolation. 


2.4 Program Logic Flow Diagram 


CONTROL PROGRAM (MAIN) 
RECEIVES THE INPUT CONTROL PARAMETERS 
INCLUDING THE BOUNDARY-LAYER 
CALCULATION OPTION BL 


BL 


BODY GEOMETRY IS 
INPUT BY CARD. 
INITIAL BODY GEOMETRY 
IS FORMULATED. 


<0 


BODY GEOMETRY WITH 
DISPLACEMENT THICKNESS 
ADD-ON IS INPUT FROM 
TAPE AND BODY GEOMETRY 


IS FORMULATED 


a {SOLVES FOR THE SOURCE DENSITY] 


COMPUTES THE COMBINED FLOWS 
AND OUTPUTS THE FINAL RESULTS 


TERMINATES 


SETS UP DATA TO BE INPUT 
TO BOUNDARY-LAYER PROGRAM 

AND CALLS THE BOUNDARY-LAYER 
PROGRAM TO GIVE 6* 


COMPUTES 3Q/a3S USES: “THE 6% 

FOR EACH ON BODY =2 PLUS THE ORIGINAL 
CONTROL POINTS BL X,y,Z TO GENERATE 
AS A NEW RIGHT = ce NEW X,Y,2 


SIDE 


RESETS BL = -BL 
RESETS BL = -BL 


2.5 Program Options 


Listed below are the options the program provides. Explanation and the 
option selector (as input variable) are given along with each option. 


1. Input moment origin 
(M@MENT ) 
Program will use an input moment origin to calculate the 
moment components in potential flow part of the program. If chosen, 
set M@MENT=1] and input the x,y,Z components of the moment origin. 


2. Input Kutta points for wake-tangency Kutta conditions 
(KUTTA) 
When this option is used, set KUTTA=1] and input the x,y,z components 
of the Kutta points as wel] as the Kutta normals. 


3. Input off-body points 
(NOFF ) 
To input off-body points, set the input parameter NOFF=1 
and input the x,y,z components of the off-body points. 


4. Partial execution — basic geometry output only 
(LIST) 
For the purpose of examining some input qeometry and not go 
throuch the whole program, set LIST=1 and the program will 
stop after body-element formation. 


oa Intermediate Print, velocity matrix 
(MPR ) 


Input MPR=1 prints the source velocity matrix Vj 
MPR=2 prints the normal velocity matrix AG j 
MPR=3 prints the dipole velocities. 


6. Intermediate print, geometry quantities and sigma matrix 


(1QUT ) 


When I@UT#0, the element geometry quantities and the solution of 


the source density matrix will be printed. 


10. 


11. 


Zs 


Lee 


Ignored element option 

(1G) 
When IG#0, input the element numbers to be ignored. The program 
will set the source velocity on these elements to zero. This 
special option is for a nonlifting section intersecting a lifting 
section. 


Special last wake element calculation 
(LASWAK) 
Input LASWAK#40 and the program will use special formulas for semi- 


infinite last wake element. In most cases, this option is selected. 


Multiple angles of attack 

(IATACK) 
If user only runs the potential flow part of the program, he may 
run it with more than one angle of attack. Set IATACK = number 
of angles of attack to be input. Maximum number allowed = 10. 
When using both potential flow and boundary-layer program, only 1 
angle of attack is currently allowed. 


Piecewise linear vorticity option 
(IWIDTH) 
If this is used, set IWIDTH#0 and input the strip widths. 


Step function vorticity option 
(ISAVE) 
Setting ISAVE#0 will use step function vorticity. 


Zero lift option 
Setting ISAVE=9 will give the zero lift output in addition to the 
normal output. This option works only when the potential flow part 
of the program is run alone. 


Symmetry features 

(SYM1, SYM2) 
Inputting SYM1#0 will treat the input as one plane of symmetry. 
With SYM140, SYM2#0 implies 2 planes of symmetry in the potential- 
flow program. 


14. Parabolic vorticity option 
(IPCV) 


If this option is used, input IPCV40. Otherwise, the program 
uses constant vorticity. 


15. Boundary-layer calculation option 
(BL) 
BL=0.0 executes only the potential flow part. 
BL=1.0 uses blowing method. 
BL=2.0 uses surface displacement method. 


16. Extra strip option 
(1XFLAG ) 
IXFLAG=0 no extra strip input. 
IXFLAG=1 treats the first input strip as an extra strip. 
IXFLAG=3 treats the last input strip as an extra strip. 
IXFLAG=2 both the first and last input strips are extra strips. 


17. Smoothing option for boundary-layer calculation 
(ISM) 
Setting ISM=] will generate 50 boundary-layer stations if the 
potential-flow input chordwise points were less than 50 on each of 
the surfaces (upper and lower). 


18. Skip the boundary-layer calculation for the nonlifting section 
(NQCAL ) 
Setting N@CAL>O will bypass the boundary-layer calculation for the 
nonlifting section. If N@CAL=0, the program assumes that the non- 
lifting section will have the boundary-layer calculation. 


2.5.1 Combination of Vorticity and Kutta Condition Option 


There are two options for selecting the Kutta condition as stated before. 
They are: 


Equal pressure Kutta condition ISAVE 
Wake tangency Kutta condition KUTTA 


These two input parameters are mutually exclusive. In other words, when 
ISAVE=0, KUTTA=1 or ISAVE=1, KUTTA=0. There are also four choices 
of vortices, namely: 
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Piecewise linear vorticity option IWIDTH=1 


Step function vorticity option IWIDTH=0 
Parabolic vorticity option IPCV=1 
Constant vorticity option IPCV=0 


If we list a table of these parameters, it can be readily seen that user may 
have eight choices of how to run his case in the potential flow part of the 
program. These choices are shown below: 


KUTTA ISAVE IWIDTH IPCV 
] 0 0 0 


o-oo 0 
=o -o-jo-+ 
—" -—- oo o—- —- oOo 
—------+=-o0o00 


2.5.2 One Point per Card Input Option 


In inputting the body points, the user has the option of punching two points 
per card or one point per card. These input formats will be explained in the 
next section. The normal way to input is two points per card. When one point 
per card is input, the user must punch the characters "QNE" on the title 
card from card column 69 to card column 72. For two points per card, this 
field is left blank. 


2.6 Program Restrictions 


2.6.1 Input Body Size 

The total number of sections input is restricted to 10, with a maxi- 
mum of 50 strips for the case. The maximum number of basic surface 
elements (including lifting and nonlifting) that can be formed from 
the input points is 1000. So if a user desires to use the full canac- 
ity of 50 strips, he may only have 20 elements on each, since 20 elements 


x 50 strips = 1000 elements. But if he chose to have 60 elements on each 
strip, then the maximum number of strips he may use is 1000 elements = 

60 elements = 16 strips. These size restrictions apply to all cases; lifting, 
nonlifting and those with both flows. 


2.4 - Onints 


- ta condition is applied by a condition of flow tangency in the 
wake (..  ,  .cing KUTTA>0), points and normal vectors must be input. 
These points are denoted Kutta points. Presently, the program allows a total 
of 50 Kutta points to be input (i.e., one for each lifting strip). The x,y,z 


coordinates of the Kutta points and the components of the Kutta point normal 


vectors are input. 


2.6.3 Off-Body Points 


When off-body points are input, the value of N@FF must be greater than 
zero. Inputting the off-body points is quite similar to that of the basic 
body points. A status flag of 3 indicates the end of the off-body point 
input. Maximum number of off-body points is 100. 


2.6.4 Angles of Attack 


User may input a total of 10 angles of attack if he uses only the 
potential-flow calculation part of the program (by inputting BL = 0.0). If 
potential-flow and boundary-layer calculations are needed, this program cur- 
rently allows only one angle of attack to be input. 


2.6.5 Nonlifting Section 


The order of input to this program should be nonlifting sections first, 
followed by the lifting sections. User may not run nonlifting sections only 
in this program, but he may run lifting sections without nonlifting sections. 
To obtain the zero lift solution, it is necessary to submit the body as a 
lifting body, then set the input parameters ISAVE=9 and BL=0; this will 
run the potential flow part alone and give the zero lift answer after the 


lifting final output. 


3.0 INPUT INSTRUCTIONS 
Some basic data inputs are necessary for the program while others are 
optional depending on the user's setup of the control card. Basically, the 


input consists of three groups. Listed below are these input groups with their 
associated card types. Each card type will be detailed in the following sections. 


3.1 Input Data Groups 


3.1.1 Title and Case Execution Information Group 


Group oe Card Type ID Remarks 
] ] Title 1 card only 
7] 2 Option control card 1 card only 
] 3 Angle of attack cards at least 1 card 
] 4 Moment origin cards at least 1 card 
(optional) 
] 5 Lifting section information cards at least 1 card 
] 6 Ignored element cards at least 1 card 
(optional ) 
] 7 Lifting strip width cards at least 1 card 


(optional) 


There are 7 types of cards, but for general use only the 4 non-optional cards 
are needed. 


3.1.2 Body Geometry Data Group 


Card 
Group Type Card Type ID Remarks 
2 ] Basic body point input cards more than 1 card 
2 2 Kutta point cards at least 1 card 
(optional) 
2 3 Kutta normal cards at least 1 card 
(optional ) 
2 4 Off-body point cards at least 1 card 


(optional) 


These are the 4 types of cards in this group. Usually only Type No. 1 is 
needed. 


3.1.3 Boundary-Layer Program Required Data Group 


Card 
Group Type Card Type ID Remarks 
3 ] Boundary-layer card 1 card only 
3 2 Transition location cards at least 2 cards 


These two types of cards are both needed. The number of transition location 
‘cards = 2 x total number of lifting strips. 


3.2 Input Card Sequence — Variables and Formats 


The input data card sequence is indicated below first with its variables 
and formats, then with a more detailed explanation of each variable. 


3.2.1 Condensed Form 


1. Title card (required) 
Format — (18A4) 
Variable — TITLE 


a: Option control card (required) 


Format — (A4, 1313, 3F4.0, F12.0) 


Variables — CASE, LIFSEC, M@MENT, KUTTA, NOFF, LIST, MPR, I@UT, IG, 
LASWAK, IATACK, IWIDTH, ISAVE, IPCV, SYM], SYM2, BL, FC 


3. Angle of attack card (required) 


Format — (6£10.0) 


Variables — ALPHAX(I), ALPHAY(I), yee 
(I goes from 1 to  IATACK 


4, Moment origin card (optional) 


Format — (6E10.0) 
Variables — MRIGNX, @RIGNY, PRIGNZ 


5. Lifting section information card (required 
Format — (1614) 


Variables — NS@RCE(J), NWAKE(J), RSTREE(?) IXFLAG (J) 
(J goes from 1 to LIFSEC 
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Vhs 


12. 


Ignored element card (optional) 


Format — (1614) 


Variables — IG1({I,J), IGN(I,J) (I goes from 1 to NSTRIP(J)) 
(J goes from 1 to LIFSEC) 


Lifting strip width card (optional 
Format — (6E10.0) 


Variables — WIDXTR(1,J), WIDTH(I,J WIDXTR(2,0) 
(I goes from 1 to atari oy) (J goes from 1 to LIFSEC) 


Basic body point input cards (required 


Either one of these two formats 
Format — (3£10.0, 211) ([2(3E10.0, 211)] 
Variables — X,Y,Z, STATUS, LABEL 


Kutta point cards (optional ) 


Format — (6E10.0) 


Variables — CUTTAX({I), CUTTAY({I), CUTTAZ(I), CUTTAX({I+1), CUTTAY(I+1), 
CUTTAZ(I+1), (I goes from 1 to KUTTA) 


Kutta normal cards (optional 
Format — (6F10.0) 


Variables — CNX(I), CNY(I), CNZ(I), CNX(I+1), CNY(I+1), CNZ(I+1) 
(I goes from 1 to KUTTA) 


Off-body point cards (optional 
Format — [2(3E10.0,11)] 


Variables — XOFF(I), YOFF(I), ZOFF(I), STATUS, X@FF(I+1), YOFF(I+1), 
ZOFF(I+1), STAT (the Tast off-body point should have 
a status = 3) 


Boundary-layer run card (required) 


Format — (2F15.0, 315) 
Variables — UI, RI, ISM, IBETA, N@CAL 


Transition location cards (required 
Format — (15) 
Variable — NTR 
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3.2.2 Detailed Explanation of the Input Data 


1. Title Card —Group 1, Type 1. Required 
Card Format — (18A4) 


Column Format Code Explanation 
1-68 17A4 TITLE 68 character description of the case to 
be run 
69-72 A4 TITLE(18) If imputting body points one potnt per 
card, punch" QNE", otherwise, Yeave 
blank. 


2. Option Control Card —Group 1, Type 2. Required 
Card Format — (A4, 1313, 3F4.0, F12.0) 


Column Format Code Explanation 

1-4 A4 CASE 4 characters used as case ID 

5-7 13 LIFSEC Total number of lifting sections input 
8-10 I3 MOMENT A nonzero number indicates coordinates 


of the moment origin will be input. 


11-13 13 KUTTA Total number of Kutta points to be input. 
If KUTTA#0, Kutta points and normals 
must be input after basic body points 
and c.c.40 on this card must be zero. 


14-16 13 NOFF Off-body point input option. A nonzero 
number indicates there will be off-body | 
points input. 


17-19 I3 LIST Case execution flag. LIST=0 implies 
full execution. LIST#0 will have 
partial execution and terminates after 
basic geometry is printed. 


20-22 13 MPR Matrix print flag 
MPR=0 no matrix print 
MPR=1 Vij matrix print 
MPR=2 Aij and sigma solution print 
MPR=3 onset flow and sigma print 


23-29 13 IQUT IQUT=0 will print all the geometric 
quantities associated with each 
element 

| 
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Column Format Code Explanation 

26-28 13 IG Ignored element option. If IG#0, card 
type #6 (ignored element card) will be 
needed. 

29-3) 13 LASWAK A nonzero number means to use the 


formulas for the special semi-infinite 
last wake element. 


32-34 13 IATACK Number of angles of attack to be input. 


35-37 be IWIDTH IWIDTH=0, step function option is used, 
and card type #7 (lifting strip width 
cards) are not needed. IWIDTH#0, piece- 
wide linear option is used, and card 
type #7 is required. 


38-40 13 ISAVE ISAVE=0, Kutta points will be used. 
ISAVE#0, equal-pressure Kutta condi- 
tion option is used and no 
Kutta points are needed. 
ISAVE=9, will also give zero lift. 


41-43 ig IPCV IPCV#0 will select parabolic vorticity 
option. 

44-47 F4.0 SYM] +1.0 or -1.0 is used for 1 plane of 
symmetry (negative sign implies anti- 
symmetric). 

48-51 F4.0 SYM2 +1.0 or -1.0 is used for two planes of 


symmetry with SYM1#0. 


52-55 F4.0 BL Boundary-layer computation option. 

BL=0 No boundary-layer calculation, 

inviscid flow only. 

BL=1.0 Use boundary-layer calculation, 
then use blowing method for the 
final computation. 

BL=2.0 Use boundary-layer calculation, 
then use the surface displace- 
ment method to compute the final 
flows. 


56-72 F12.0 FG Foot converter to convert the input 
coordinates to the units of feet. 


3. Angle of Attack Card — Group 1, Type 3. Required 
Card Format — (6£10.0) ix 


Column Format Code 

1-10 E10.0 ALPHAX 
11-20 E10.0 ALPHAY 
21-30 E10.0 ALPHAZ 


Explanation 


x-component of the first uniform onset 
flow. 


y-component of the first uniform onset 
flow. 


z-component of the first uniform onset 
flow. 


For more than J] uniform onset flow, 

repeat the same format from c.c. 31-60. 
Two uniform onset flows (angles of attack) 
per card. May use as many cards as 

needed to satisfy JATACK. 


4, Moment Origin Cards — Group 1, Type 4. Optional 


Card Format — (6£10.0) 


Column Format Code 

1-10 E10.0 QRIGNX 
11-20 E10.0 @RIGNY 
21-30 E10.0 MRIGNZ 


5. Lifting Section Information Card 


Card Format — (1614) 


Column Format Code 
1-4 14 NSORCE (J) 
5-8 14 NWAKE(J ) 
9-12 14 NSTRIP(J) 
13-16 14 IXFLAG (J) 


Explanation 


x-coordinate of the moment origin. 
y-coordinate of the moment origin. 
Zz-coordinate of the moment origin. 


-— Group 1, Type 5. Required 


Explanation 


Number of on-body elements (including 
ignore-elements) on each lifting strip 
in the lifting section J. 


Number of wake elements on each lift- 
ing strip in the lifting section J. 


Number of lifting strips (including 
extra strip) in the lifting section J. 


Extra strip option. 

IXFLAG=0 no extra strip input 

IXFLAG=1 first input lifting strip of 
section J is an extra strip 

IXFLAG=3 last input lifting strip of 
section J is an extra strip 

IXFLAG=2 implies both first and last 
strips of a lifting section 
J are extra strips 
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Column Format Code Explanation 
Repeat the above format again from c.c. 
17 to 64. Four sections per card. 
Subscript J goes from 1 to the total 
number of lifting sections (LIFSEC). 
Ignored Element Card — Group 1, Type 6. Optional 


Card Format — (1614) 


Column Format Code Explanation 
1-4 14 IG1(I,J) Punch a number K — it means the Kth 


element in a lifting strip that is used 
as the beginning of the ignored element 
group. 


5-8 14 IGN(I,J) Punch a number M — it means the Mth 
element in a lifting strip that the 
ignored element group ends. For example, 
assume a lifting strip has 10 source 
elements, and the 3rd, 4th, 5th elements 
are the ignored elements. In this case, 
IGl=3, IGN=5. I = strip indix, J = sec- 
tion index. Repeat the same format on 
the same card for the next strip. 

7 strips of ignored elements per card. 
Start a new card for a new lifting 
section. 


Lifting Strip Width Card — Group 1, Type 7. Optional 
Card Format — (6£10.0) 


Column Format Code Explanation 
1-10 E10.0 WIDXTR(1,J) Width for the first extra strip of a 


lifting section. Leave blank if there 
is no first extra strip. 


11-20 E10.0 WIDTH(I,J) Width value for the first lifting strip. 
Repeat the same format for successive 
strips until all lifting strip widths 
are punched. If there are more than 
6 values, use another card. I = strip 
index, J = section index. 


21-30 E10.0 W'DXTR(2,J) The last width value is the last extra 
strip width. If no‘last extra strip, 
leave this field blank. Start a new 
card for a new lifting section. 
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Basic Body Point Input Cards — Group 2, Type 1. Required 
Card Format — [2(3E10.0,211)] if TITLE(18) jis blank. 
Card Format — (3E10.0,211) if TITLE(18) is " QNE". 


Column Format Code Explanation 
1-10 E10.0 X X,y,Z coordinates of a point on the 
11-20 E10.0 Y line to be read in to form an element 
21-30 E10.0 Z 
31 1] STATUS STATUS=2 new section 


STATUS=]_ new line 
STATUS=0 same line 
STATUS=3 end of all body point input 


32 a LABEL LABEL=1 lifting section 
LABEL=0 nonlifting section 
ae a X,y,Z coordinates of the next point. 
53-62 F10.0 71 Use only if 2 points per card option 
‘ was selected. Status & label have 
63 1 STATUS same meanings as above. 
64 1] LABEL 


Kutta Point Cards — Group 2, Type 2. Optional. 
Card Format — (6£10.0) 


Column Format Code Explanation 

Bg oe palit X,y,Z coordinates of the Kutta point 
21-30 E10.0 cuTTAz(r) controls 

ie ane et a ee X,y»z coordinates of the next Kutta 
51-60 E10.0  CUTTAZ(I+1) POTts 


Two Kutta points per card. Repeat same 
card until al] Kutta points are input. 


Kutta Normal Cards — Group 2, Type 3. Optional. 
Card Format — (6E10.0) 


Column = Format. Code Explanation 
1220 ae in X,¥,sZ components of the Kutta point 
21-30 E10.0 CNZ(1) normal vector 


2] 


Column Format Code Explanation 


ae A etelias X,y,Z components of the next Kutta 
51-60 10.0 CNZ(I+1) point normal vector 


Repeat the same format for the next 
card if needed, until all Kutta point 
normals are input. 


Off-Body Point Cards — Group 2, Type 4. Optional. 
Card Format — [2(3E10.0,11)] 


Column Format Code Explanation 
1-10 E10.0 XMFF (I) ; 
11-20 F10.0 Y@FF (I) ae components of the input off-body 
21-30 E10.0 ZOFF (1) P 
3] 1] STATUS Punch 3 for the last off-body point, 
otherwise leave blank 
32-41 E10.0 X@FF(I+1) 
42-5] E10.0 YOFF(I+1) Same meaning as above 
52-61 E10.0 ZOFF (I+1) 
62 1] 


Repeat the format of this card as many 
time as needed for all the off-body 
points. The last card must have a 3 
punched in either column 31 or 62 
depending on where the point ends. 


Boundary Layer Run Card — Group 3, Type 1. Required. 
Card Format — (2F15.0,315) 


Column Format. Code Explanation 
1-15 F15.0 UI Reference velocity, (ft/sec) 
16-30 F15.0 RI Reynold's number, (ft!) 
31-35 LS ISM ISM40 will generate more boundary-layer 


calculation stations from the potential - 
flow output stations. 


36-40 15 IBETA IBETA#0 uses alternate way to compute 
a. Usually this flag should be set to 0. 


41-45 15 N@CAL NOCAL#0 will suppress the boundary-layer 
calculation for the nonlifting section. 
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13. Transition Location Card — Group 
Card Format — (15) 
Column Format Code 


1-5 Ie NTR 


3.3 Input Card Deck Setup 


The following illustration shows 


TITLE CARD (1,1) 


3, Type 2. Required. 


Explanation 


Transition location of the input bound- 
ary layer strip. If this value is 
unknown, user may punch a large value 
and the program will calculate the 
transition location. Two transition 
location cards per lifting strip, one 
for the lower surface and one for 

the upper surface, are input in that 
order. The number of transition loca- 
tion cards needed = 2 x the total number 
of lifting strips. 


the arrangement of the input data cards. 


TRANSITION LOCATION (3,2 
BOUNDARY -LAYER_RUN CARD 


First index on the upper right corner indicates the Group Number, second one 


is the Card Type Number. 
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4.0 INPUT LOAD SHEETS 


The following input data sheets are designed to assist users loading 
their data to be key punched. The Group No., Card Type No., and Sequence No. 
are put in for reference only and the user need not fill them in. 
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5.0 OUTPUT EXPLANATION 


The output of each case is composed of four parts. Part one is the 
initial potential flow run. Part two is the output from the simulation program. 
Part three has all of the boundary-layer-calculation output. Finally, part 
four is the final output of the potential-flow calculations. 


In addition to the above normal output, there are some optional inter- 
mediate outputs, which will be explained after the normal output. 


5.1 Results of the Initial Potential-Flow Calculations 


5.1.1 Basic Data 


Symbol Definition 

N An integer index to identify the strip number of the section. 
M An integer index to identify the element number of the strip. 
X x-coordinate of the element's corner point. 

Y y-coordinate of the element's corner point. 

Z z-coordinate of the element's corner point. 

NX X-component of the unit normal vector. 

NY y-component of the unit normal vector. 

NZ z-component of the unit normal vector. 

X0 x-coordinate of the centroid of the element. 

YO y-coordinate of the centroid of the element. 

Z0 z-coordinate of the centroid of the element. 

D The distance the four corner points are projected to become coplanar. 
T Maximum diagonal length of the projected element. 

A Area of the element. 

NLIF Nonlifting section element. 

LIFT Lifting section element. 


XTRA Extra strip element. 


Symbo1 Definition 


WAKE Wake element. 


Note: This will be the only output if "LIST" is input with a nonzero value 
on card #2, 


5.1.2 On-Body Solutions 
Symbo] Definition 


N An integer index to identify the strip number of the section. 
M An integer index to identify the element number on the strip. 
XO x-coordinate of the centroid of the element. 
YO y-coordinate of the centroid of the element. 
Z0 z-coordinate of the centroid of the element. 
VX x-component of the total flow velocity. 
VY y-component of the total flow velocity. 
VZ z-component of the total flow velocity. 
VT Magnitude of the total velocity. 
Magnitude of the total velocity squared. 
Pressure coefficient, equals 1 —VTSQ. 
x-component of the direction cosine of the velocity vector. 
y-component of the direction cosine of the velocity vector. 
z-component the direction cosine of the velocity vector. 
x-component the unit normal vector. 
y-component the unit normal vector. 
z-component the unit normal vector. 
Surface source density. 
Dot product of the velocity and the unit normal. 


Area of the element. 


Note: 


ra 


Symbol] 


KUTTA 
POINTS 


XO 
YO 
Z0 
VX 
VY 
VZ 
VT 
VN 
cP 
DCX 
DCY 
DCZ 
NX 
NY 
NZ 


5.1.4 


Symbo] 
POINTS 


XO 


1. The components listed for the forces and the moments are in the 
Xyy,Z order. 


2. If there are ignored elements input, the solutions for the ignored 


elements will not 
in this case will 


be printed. The output of the on-body solutions 
show less elements than other strips which do not 


contain the ignored elements. 


Kutta Point Solutions 


Definition 


Identifies the Kutta point number. 


x-coordinate of the Kutta point. 


y-coordinate of the Kutta point. 


z-coordinate of the Kutta point. 


x-component of the 
y-component of the 


z-component of the 


total velocity. 
total velocity. 


total velocity. 


Magnitude of the total velocity. 


Dot product of the 


velocity and the unit normal. 


Pressure coefficient. 


X-component of the 
y-component of the 
z-component of the 
X-component of the 
y-component of the 


z-component of the 


direction cosine of the velocity vector. 
direction cosine of the velocity vector. 
direction cosine of the velocity vector. 
Kutta normal vector. 
Kutta normal vector. 


Kutta normal vector. 


Off-Body Point Solutions 


Identifies the off- 


Definition 


body point number. 


x-coordinate of the off-body point. 
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Symbol Definition 
YO y-coordinate of the off~body point 

Z0 z-coordinate of the off-body point 

VX x-component of the total velocity 

VY v-component of the total velocity 

VZ z-component of the total velocity 

\T Magnitude of the total velocity 

VTSQ Magnitude of the total velocity squared 
cP Pressure coefficient 

DCX X-component of the direction cosine 
vCY y-component of the direction cosine 
DCZ z-component of the direction cosine 


5.1.5 Bound Vorticity Coefficient Solution 

Symbol Definition 

STRIP NQ. Identifies the order of the strips in the lifting section 

B(STRIP) The bound vorticity coefficient (chordwise dipole derivative) 
associated with each strip. 

5.1.6 Intermediate Print 


This section explains the various messages, notes and quantities that are 
printed, other than the final output explained before. 


We Control Flags Print 


This is the first page printout of all the input control flags and 
all the uniform onset flows. This part of the output enables user 
to check if the control variables are input properly. 


2. Time Print at Various Steps During Execution 


When the computer run starts, time is initialized to zero. At various 
parts of the program during execution of the run, current CPU time 

is called and printed. Subtraction of the time between steps will 
indicate the amount of CPU time used. 
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errr ener ens enri tanner ertresecnnceemeccacmaeall 


BEGIN THE SUBR@UTINE INPUT. TIME IS xxx.xxx SEC. 
Note: The time at the beginning of the input routine. 


BEGIN THE F@RMATIQN QF ELEMENTS. TIME IS xxx.xxx SEC. 
Note: The time at the start of element formation. 


FINISH THE FORMATION OF ELEMENTS. TIME IS xxx.xxx SEC. 

Note: The time at the completion of element formation including 
the printout of the basic data. 
Subtraction of the previous step time from this time should 
give the time it takes to form the body elements. 


END QF SUBR@UTINE INPUT. TIME IS xxx.xxx SEC. 
Note: The time at the end of the input routine. 


BEGINNING THE VF@RM R@UTINE. TIME IS xxx.xxx SEC. 
Note: The time at the beginning of the velocity formation routine. 


LIFTING 
NONLIFTING 


Note: This is the time required for forming either the nonlifting or 


TIME F@R THE FORMATIPN OF THE 


lifting velocity matrix in minutes. 


END OF THE VFQRM ROUTINE. TIME IS xxx.xx SEC. 
Note: This is the time at which the formation of the velocity matrix 
for the whole body is completed. 


END QF THE AF@RM ROUTINE. TIME IS xxx.xx SEC. 
Note: This is the time at which the formation of the Aij matrix 
is completed. 


THE NX N MATRIX WITH m RIGHT SIDES WAS SOLVED DIRECTLY IN xx.xxx 
MIN. 
Note: This gives the total time the matrix solution routine requires. 


Source Velocity Matrix Print 


COLUMN 
CNTRL PT 
VXS 

VYS 

VZS 
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VEL@CITY MATRIX = xx.xx MIN. 


ee, 


Explanation: The numbers printed below these headers are the source 


velocity matrix column number, the control point number, the x,y and 


z components of the source velocity matrix. 


Note: This printout will only appear when MPR=1 is input. 


Dipole Velocity Matrix Print 
STRIP 
CONTROL PT 
VXF, VXS 
V¥F, ¥¥S 
VZF, VZS 


Explanation: The numbers printed below these headers are the strip 
count (lifting), control point number, the x,y,z components of the 


first and second dipole velocity. 


Note: This printout will appear only if MPR is set equal to 2. 


Onset Flow Matrix Print 


QNSET FLOW NO 
CONTROL PAINTS 
X-FLQW 
Y-FLQW 
Z-FLOW 


Explanation: The numbers printed below these headers are the onset 
flow number, the control point number, and the x,y,zZ components of the 


onset flows. 


is input. 


Note: This printout will appear only if MPR=2 


The Dot Product (Aij) Matrix Print 
COLUMN 

Ald 
FLOW NO. 
RIJ 


Explanation: Under these headers are the column count of the AlJ matrix, 
the dot product value, the onset flow numbers and the righthand side matrix. 


Note: This printout will appear if MPR is set greater than or equal to 2. 
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7. The Sigma Matrix Print 


SOLUTI@N @BTAINED AFTER C@LS@L 
FLQ@W N®. — Indicates the onset flow number. Following it are the 
sigma solutions. 8 numbers to a line. 


Note: This matrix will be printed if MPR is not less than 2. 


8. Other Intermediate Print 

TABLE QF INPUT INFORMATION 

INPUT SECTIQ@N N@. — Section number count 

SECTI@N TYPE — Type of input section 
0 nonlifting 
1 lifting 

TOTAL N@. OF ELEMENT IN EACH SECTION 
Summarize the total element count (including wake 
element less extra strip) 

EXTRA STRIPS —0 no extra strip 
) first strip is extra 
3 last strip is extra 
2 both first and last strip are extra strips 

STRIP NO. — Indicates the strip count 

SPURCE ELEMENTS IN THE STRIP — Prints number of on-body elements in 
each strip (including ignored elements). 

WAKE ELEMENTS IN THE STRIP ~ Indicates the wake elements in the 
strip. 

TOTAL N®. OF ELEMENTS INPUT — Sums up all the elements input (not 
counting the extra strip elements). 

BODY SECTIG@N N@. = n — Indicates that the velocity formation is 
being processed for the input section n. 

TYPE = i — Indicates this section n is of the type i(i = 0 for 
nonlifting, i= 1 for lifting). 

TOTAL NO. BF PPINTS = mm — Indicates the total elements in the 
section n. 

NQ. OF STRIPS = j — Indicates the total number of strips in the 
section. 

TOTAL NO. OF CONTROL PHINTS (INCL. QOFF-BODY PPINTS) = nn — Indicates 
the total number of control points in the section. 
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LIFTING STRIP N@. j, NO. OF IGNORE ELEMENTS k — Indicates the 
lifting strip count and its number of ignored elements. 

TOTAL NO. MF ELEMENTS IN THE LIFTING SECTI@N — Indicates the total 
number of lifting elements in the lifting section. 

N@. QF FAR ELEMENTS = xx — Far field element count. 

NO. QF INTERMEDIATE ELEMENTS = xx — Intermediate field element count. 

N®. OF NEAR FIELD ELEMENTS = xx — Near field element count. 


9. Identifying the final output for certain uniform onset flow 
FINAL QUTPUT FQR THE FOLLOWING UNIFQMRM PNSET FLOW 
(EX. XXKANK, EX AXKKAK, FX. KKXK) 


These are the components of the input uniform onset flow. 


10. Identifying the zero-lift output for certain uniform onset flow. 
ZER@ LIFT QUTPUT FOR THE FOLL@WING UNIFQ@RM QNSET FLOW 
CF MIOUKK 5 MAK, $KKKKK) 
These are the components of the uniform onset flow from which 
the zero lift is calculated. 


5.2 The Simulation Program 


There are two different types of simulation program output: one is for 
blowing method, and the other is the result of the surface displacement method. 
Either type of output is printed depending upon user's choice. 


5.2.1 Before Boundary-Layer Program 


This part of the output is the same for both blowing and surface 


displacement methods. 
Symbol] Definition 


SECTION NQ. Section number of the input body 


TYPE Nonlifting section is type 0. 
Lifting section is of type 1. 


TOTAL STRIPS Total number of strips within the section. 


STRIP NQ. Strip number of the printed section 


STAGNAT IQN Element number of the strip where stagnation point is 
POINT calculated 


Symbol 


XO 
YO 
Z0 


BL. STRIP 


C. POINT 


K 
PCH@RD 
XK(K) 
ZK(K) 
S(K) 
V(K) 


Definition 


The x,y.Z coordinates of the computed stagnation point. 


Boundary-layer strip number. A value of 1 denotes it is the 
lower surface of the original strip used in the potential flow 
calculation, while a value of 2 is for the upper surface. 
Original potential-flow calculation control point numbers, 
except number one which is the computed stagnation point. 


The x,y,Z components of the control point. 


Surface velocity at the control point. 
Arc length between consecutive control points. 
Boundary-layer station numbers if smoothing option is chosen, 


The percent chord associated with the smoothing stations. 
Note: This is a prescribed array in the program, 


The interpolated x,z values at the smoothed boundary-layer 
stations. 


The interpolated arc lengths between successive stations. 


The interpolated surface velocities. 


5.2.2 After Boundary-Layer Program 


Symbol 
I 


SM(1) 
VB(I) 


DEL (1) 


Definition 
Potential-flow calculation station number 


Back interpolated arc length values from the smoothed boundary- 
layer stations. 


Back interpolated surface velocities from the smoothed boundary- 
layer stations. 


Back interpolated 6* values from the smoothed boundary-layer 
stations. 
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5.2.3 Surface Blowing Simulation Output 


1. THE NEW RIGHT-HAND SIDE 
These are the generated blowing normal velocities to be used in 


solving for the final potential-flow solution. These velocities 

for each strip are printed after the strip has finished its boundary- 
layer calculation, and, finally, a whole set of blowing normal 
velocities is printed before the matrix solution routine is called. 
This set of blowing normal velocities will be the right-side colum 
matrix [C] in the matrix equation [A] cy"! = [B], where 


[A] is the [Aij] matrix from the initial potential-flow 
calculation, 


[B] is the source density column matrix to be solved. 


a S@LUTIQ@N MPBTAINED AFTER COLSO@L 
FLOW NO. 1 
This is the solution matrix [B] (the source density matrix) obtained 
after the matrix solution routine is called. 


3. CHECKING THE SIGMAS AFTER B. LAYER LINK 
SIGMAS QF THE UNIF@RM FLOW NO. 1 
First the original o solution of the uniform flow is printed. 
This enables the user to check the same uniform flow solution (from 
initial potential-flow calculations) as brought back. Then the 
second part with the same message printout is the uniform flow oa 
solution plus the solution obtained from 2. It is the final o 


result to be used in computing the final flow. 


4. — B-ARRAY 
The new bound vorticity coefficient associated with each strip. 


5. Boundary-layer calculation summaries for each lifting strip. 
A table, which consists of the input to the boundary-layer program 
and the output 45*'s for each lifting strip, is provided. Explana- 
tions of these parameters are given below. 
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symbol Definition 

STRIP NQ. Lifting strip number. 

C.PQINT Control point number in the strip 

: X,Z coordinates of the control point 

Vv Surface velocity at the control point 

S Arc length between two control points 

DEL Displacement thickness 6* at control points 

DOQDS 3q/3S computed for each control point, where q = 6*V. 


When 3aq/3S over the entire body is computed, it is 
used as a new right-hand side in the original matrix 
equation to be solved for density distribution. 


5.2.4 Surface Displacement Simulation Output 


After both the lower and upper surfaces are computed in the boundary-layer 
program, a summary of the input and output variables are printed in one table. 
The following list shows the symbols and their definitions. 


Symbol 


STRIP NO. 


C.PQINT 


X 
7: 


DEL 


Definition 
Lifting strip number 


Control point number in the strip 
X,Z coordinates of the control point 


Surface velocity at the control point 
Arc length between two control points 


Displacement thickness 6* at control points 


Another output is a table of §*. This output is solely for the purpose 


of checking the 46*'s on the external unit before they are used in generating 


new coordinates. 


Symbo] 
DEL STAR 


K 


KI 


Displacement thickness 4* 
A do-loop counter counting the number of strips 


A counter that was stored on the tape to identify the tape 
number. Kl should always be equal to K. 
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5.3 Boundary-Layer Program Output 


The output of this program consists of three parts: The input body 
geometry data, the computed station data, and the output summary at each 
station. 


Listed below are the output variables and their definitions in the order 
of their appearances in the output. 


Symbol] Definition 
TRFLAG Transition location flag. The value of this flag is set equal 


to 1. The program will calculate the transition location 
if the user inputs a number larger than the number of 
boundary-layer calculation stations for the value of NTR. 


TRINT Transition control flag. This flag is set equal to zero which 
means the transition is instantaneous. 

TVC = 0, not applicable 

SH@RTP Print control flag. This flag is set equal to 1 in the 


program. This means the program will use short printout (no 
velocity or enthalpy profiles printed). 


X/C Nondimensionalized x-coordinates. 

Y/C Nondimensionalized y-coordinates (it is actually the z-values 
from potential-flow program) 

X x-coordinate in feet 

S Surface distance in feet 

S/C Nondimensionalized surface distance 

H] First increment of » calculation 

C Chord length in feet 

PRO Laminar Prandtl] number 

K Variable grid factor. It is set equal to 1.14. 

RH@REF Freestream reference density, slugs/ft> 

SWEEP Sweep angle. Not used in this version. 

KK Flow index KK = 0.0 for 2-D flows 

MUREF Freestream dynamic viscosity uu, Ibe sec/fté 
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Symbol 


HE 
EPS] 
VREF 
TREF 
REY 
MREF 
RO/C 
TW 
UE 
PE 
FW 


ALPHA] 
QW 


cP 
MUE 


FPW 
BETA 
ALPHA2 
RR 

ME 

TE 


SQUIG 


THETA 


Definition 
Total enthalpy at the edge of boundary layer, (ft/sec)*. 
Converge criteria for laminar flow. Set equal to 0.005. 
Reference velocity, ft/sec. 
Reference temperature T , °R. 
Reynolds number based on reference conditions. 
Freestream Mach number . 
It is the z values from the potential-flow program. 
Not applicable here. 


u velocity at edge of boundary layer, ft/sec. 


e? 
Pe, pressure at edge of boundary layer, Ib,/ft®. 


Ts the transformed stream function at the wall 


1 / Pwr 
ty we oF paca, 
" (2c) Petette 


Not applicable to this version 


Pressure coefficient 


Local dynamic viscosity jie, at edge of boundary layer 
Ibe + sec/ft? 


f! 
W 


8 = (2€/u,)(du,/de) 
Not applicable in this version 


Transformed x-coordinate ¢& 
x 


& s 
, S ere'e dx 


0 


Momentum thickness, 98, ft 
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Definition 


Boundary-layer displacement thickness 


Value of the local skin friction coefficient, Cr 


Ce = ty/ (1/20 4u6) 


ig 

Total enthalpy ratio H/H, 

Number of points taken through the boundary layer. 
RX Reynold's number based on x 
RTHETA Momentum thickness Reynolds number, R, 
H Boundary-layer shape factor, H = 6*/6 
CFA Integrated skin-friction coefficient, Cr 

qy? slope of gy asa function of n at the wall 


ST Stanton number, not used in this version 


ETAINF No 
5.4 Final Output 

The final output format of the program is identical to the initial 
potential-flow calculation output. When the blowing method is used the 
original coordinates will be output again with the new potential flow solu- 
tion. For the surface displacement method, the new body coordinates are 
printed with the final potential-flow solution. 


6.0 ERROR MESSAGES 


6.1 Inviscid Flow Program 


a. Message: MISMATCH @F ELEMENTS IN A LIFTING STRIP IS DETECTED. 
ELEMENTS F@RM = xxx, ELEMENTS INPUT = xxx, CQMPUTATION 
TERMINATED. 


Cause of error: Inconsistent input data. The program sums up the 
number of on-body elements plus the wake elements. This 
sum does not match that of the input elements formed from 
the input coordinates. 


Action: Check the lifting body information card and the input basic 
body points coordinate cards. 
Number of points input should equal the number of elements 
plus 1 on each strip. 


For example. If in a lifting section, each lifting strip 
consists of 10 on-body elements and 1 wake element, the 
total element number is 11, and there should be 12 points 
input. Each point has x,y,z coordinates. 


6. Message: ERR@R IN IGNORED ELEMENT CQUNT xxx, SH@ULD BE xxx. 
Cause of error: Erroneous input of the ignored element information. 


Action: Check the input data on group 1, type 6 cards. Make sure 
the ignored element information is properly input. 


c. Message: LABEL ERROR IN NONLIFTING VFQRM. 


Cause of error: Reading in wrong part of data from unit 4. The 

geometric quantities of each element are stored in unit #4. 
These data must be read in during the velocity formation 
routine. In front of these geometric quantities is a label 
which indicates they are lifting section data or nonlifting 
data. If computation is in the nonlifting section and the 
system is reading in lifting section data, the above mes- 
sage will be displayed. 
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Action: Check that the number of lifting strips input 
is the actual number of strips input. 


Message: ERROR IN VF@RM. THE ELEMENTS FQRMED DM NOT CMRRESPOND TP 
THE NO. OF BODY ELEMENTS. 


Cause of error: The input number of elements does not correspond 
to the actual number of elements being formed. 


Action: Check lifting body information inout, and also the non- 


lifting body points, if any. 


Message: AFTER xxx ITERATI@NS, DELTA B STILL DID NOT CONVERGE TQ 
THE GIVEN CRITERI@N/LARGEST DELTA B = + x.xxxxxxE*xx/ 
PROGRAM PROCEEDS WITH THE MOST CURRENT V@RTEX STRENGTH. 


Cause of error: In the stepwise routine, after 100 iterations, the 
AB value is still bigger than the given criterion, so the 
max. B value is used for further calculation. (B = Bound 
vorticity coefficient. ) 


Action: Check the basic body coordinate input. 


Message: XXX QNBADY POINTS MISSED, EXECUTI@N TERMINATED. 


Cause of error: In the final output, the on-body element count 
is incorrect. 


Action: Check the number of source elements input. 


Message: XXX KUTTA P@INTS MISSED, EXECUTION TERMINATED. 


Cause of error: In the final output, the Kutta point count is 
incorrect. 


Action: Check the total number of Kutta points input. 


Message: XXX OFF-BADY PAINTS MISSED, EXECUTI@N TERMINATED. 


Cause of error: In the final output, the off-body point count is 
incorrect. 


Action: Check the number of off-body points input. 
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ij. Message: THE MATRIX OF XXXXX ELEMENTS EXCEEDS xxxxxx. 


Cause of error: The size of the matrix is too big for the work area 
prepared for the CQ@LSPL routine. 


Action: Change ISIZE in the main program. 


6.2 Boundary-Layer Program 
a. Message: ***ITERATIQNS EXCEED THE ALL@WABLE LIMIT*** 


Cause of error: In subroutine BOUNDL, if the number of iterations 
exceeds 9, this message will be printed and job stops. 


Action: None. 


b. Message: NEGATIVE VEL@CITY IS FQUND, ITERATIQN ENDS 


Cause of error: During iteration at a given station, a negative 
velocity is found, and run stops. This error message is 
printed in subroutine B@UNDL. 


Action: None. 


c. Message: **ERROR — NO INPUT FOR EITHER VREF QR MREF 


Cause of error: Both the reference velocity UI and the reference 
freestream Mach number are missing from the input. 


Action: Check the input of VREF, make sure it is a nonzero number. 


d. Message: **ERR@R — XI AT STATION 1 NE QR GT O. 


Cause of error: The first point of the surface distance must be 
equal to 0. 


Action: Check the input S array. If S(1)#0. Change it and 
rerun the job. 
e. Message: **ERR@R — FP PROFILE IS NEGATIVE AT I = xxx 


Cause of error: At the printed station number, the f' profile is 
found to be negative. 


Action: None. 


Message: **ERR@R — IMAX(xxx) EXCEEDS 100 — IMAX = xxx. 


Cause of error: The number of iteration at the printed station 
number exceeds 100. 


Action: None. 


Message: ******#*k*CASE TERMINATED****** 4% 4% 


Cause of error: Any occurrence of the error messages in the 
boundary-layer program. 


6.3 Simulation Programs 


a. 


Message: N@ SEPARATIQN CAN BE F@UND F@R THIS STRIP 


Cause of error: Using the potential-flow calculation results, the 
simulation program can not locate the stagnation point to 
separate the lifting strip into two boundary-layer strips. 


Action: Check the input and output data of the potential-flow 
program. 


Message: TOTAL PAINT IN SETUP NOT EQUAL TQ TOTAL CONTROL POINT, 
PROGRAM ENDS FOR C@RRECTION 

Cause of error: In generating the new onset flow, the program checks 
the total number of control points with the control ele- 
ments that are currently being worked on. If these two 
numbers do not agree, this messaqe is printed. 


Action: Check that the total number of control elements is passed 
correctly into here, and also check to see if some of the 
control points have been skipped. 


Message: ELEMENTS INPUT T@ B.L. PR@GRAM = xxxxx. ELEMENTS RETURN 
T@ CALLING PR@GRAM = xxxxx. PR@GRAM STOPS. 


Cause of error: Inconsistent source element count. 


Action: Check that the number of source elements per strip stored 
on the external unit is agreeable with the number already 


stored in core. 


Message: STRIP CPUNT ERR@R IN UNIT 3 


Cause of error: The strip number stored on unit 3 with the 6*s is 
incorrect. 


Action: Check that the total number of lifting strips is correctly 
passed. 
Message: STRIP C@UNT DIFFERS FROM TAPE IDEL 


Cause of error: The strip counter recorded on the unit IDEL for the 
lifting section is not the same as the strip number being 
used at the time the error occurs. 


Action: Check the contents of the unit IDEL. 


Message: ELEMENT N@. DIFFERS FROM TAPE IDEL 


Cause of error: The number of on-body elements recorded on unit 3 
is incorrect. 


Action: Check the contents of the unit 3. 


Message: FINAL STRIP DIFFERS FROM TAPE IDEL 


Cause of error: The total number of lifting strip count is not 
agreeable with the strip index on unit 3. 


Action: Check the contents of unit 3. 


Message: WAKE ELEMENT C@UNT ERR@R. Kl = XXX MI = XXX EXECUTION 
STOPS WITH ERR@R CPDE 20 


Cause of error: Generated total number of on-body points plus wake 
points is not the same as stored from the initial potential 
flow calculation. 


Action: Check that the correct number of on-body elements and wakes 
are stored and passed to this part of the program correctly. 
Message: IWAKE = XXXXX IWK = XXXXX ST@P 19 


Cause of error: The number of generated wake elements is not the 
same as stored from the initial potential-flow calculation. 


Action: Check the number of elements and wakes generated in this 
part of the program. 
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7.0 DECK SETUP 


This section illustrates how the input deck should be set up for various 
types of computer runs. The overlay arrangements and the job control cards 
illustrated are good for the IBM 370/165 systems only. Proper conversion 
should be made when other systems are used. 


7.1 Overlay Structure 


Because of the size and the structure of the program, it is beneficial 
to use the overlay technique to minimize the case region when executing it. 
The overlay structure is shown with segment numbers on the following page. 
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This page shows the overlay control cards corresponding to the overlay structure 
shown on the previous page. 


I 2 3 4 $ 6 
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7.2 The Data Definition Cards (DD Cards) for the External Units 


There are 17 external units required by this program. Special attention 
should be paid to the SPACE parameter. The number of tracks denoted in the 
SPACE parameter depends on the total number of input elements that defines 
the basic body. For instance, the unit that stored the Aij matrix for a 
300-element case may require 40 tracks on a 3330 disk, the same unit when 
used to store a 1000 element case will require 320 tracks. User should esti- 
mate the track sizes before submitting a computer run. 


Example No. 1 shows the DD Cards for cases with less than 400 elements 
and example No. 2 shows them for cases of 1000 elements or less. 


7.3 Job Control Language Cards 


The IBM/370 JCL card setup for two types of computer run is shown below. 
The first example shows a COMPILE, LINKEDIT and G™ job. The second one shows 
the G®@ step alone, i.e., executing the program which is already stored on a 
disk pack. 


In these two examples, the data set name DAC.EOH3,DPM.J1HB.V2 is used. 
The volume on which the data set is stored is called CSLB30. The DD cards 
from examples 1 or 2 may be used or redefined by the user depending on the 
case size. 


7.3.1 JCL Setup for Compile, Linkedit and Go. (First line is the card column 
identifier, not a JCL card.) The deck setup for this kind of job is shown in 


example No. 3. 


7.3.2 JCL Setup for Executing the Program Already Stored on Disk. (First 
line is the card column identifier, not a JCL card.) This type of run is 


shown in example No. 4. 
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7.4 Run Time Estimation 


The estimation of run time depends a great deal on the total number of 
elements that depicts the input body, and the number of lifting strips. In 
general, the surface displacement method will run longer than the blowing method. 


The figure below shows the number of elements versus time in minutes. 
The points shown on the figure are the actual data from various runs made on 
the IBM/370 systems. Users may apply these curves as a guide for their CPU 
time estimation. 
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8.0 TEST CASE 


8.1 Test Case Input Data 


The test case is a four strip swept wing with a symmetric airfoil section 
run at 


¥, 
as 
3. 


8. 
oe 


The lifting body input consists of 4 lifting strips, and on each strip there 
are 30 on-body elements and 1 wake element. 


a = 8.22°. Listed below are the options used: 


One angle of attack 

One plane of symmetry 

One lifting section 

Special formula for last wake 
Stepwise vorticity 

Equal pressure Kutta conditions 
Case 1, blowing method 

Case 2, displacement method 


Input foot-converter 


following page. 


A complete set of data for test case No. 1 (blowina) is included in the 
following pages. 


TATACK=1 
SYM1=1.0 
LIFSEC=1 
LASWAK=1 
IWIDTH=0 
ISAVE=1 
BL=1.0 
BL=2.0 
FC=12.0 


A sketch of this body geometry 
showing the wing planform with the strip locations, and the airfoil configura- 
tion with its defining elements (without wake element) is presented on the 


For test case No. 2 (surface displacement), only the first 
two cards (group 1, type 1 and 2) from case No. 1 need to be changed, hence, 
only these two cards are listed. 


ACCOUNTED FOR 
BY SYMMETRY 


WING PLANFORM 


AIRFOIL SECTION 


Sketch of Test Case Geometry. 
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8.2 Test Case Output 
The test case output in the following pages consists of: 
a. Complete potential-flow-calculation output. 
b. One strip of boundary-layer-calculation output. 
c. Complete sets of final outputs for both blowing and surface displace- 


ment methods. 


The initial potential-flow-computation output is identical in both blowing and 


surface displacement methods, and one set of this output is included here. 


Due to the amount of printout generated in the boundary-layer calculation, 
only 1 strip of boundary-laver output is shown as an example since the other 
strips have the same output formats. 


8.2.1 Initial Potential-Flow-Calculation Output. (Same for both blowing and 


displacement methods. ) 
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8.2.2 Boundary-Layer-Calculation Output 


Strip No. 4 is shown. Also included is the simulation program output. 
This part of the output is the same for both blowing and displacment methods. 
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